FIELD OF THE INVENTION 

The object of the present invention is a method for gradually deforming an initial 
distribution of objects of geologic nature, formed by simulation of a stochastic model of 
object type, from measurements or observations, so as to best adapt it to imposed 
physical constraints of hydrodynamic type for example. 

The method according to the invention finds applications in the sphere of 
underground zones modelling intended to generate representations showing how a 
certain physical quantity is distributed in a subsoil zone (permeability for example) and 
which are best compatible with observed or measured data : geologic data, seismic 
records, measurements obtained in wells, notably measurements of the variations with 
time of the pressure and of the flow rate of fluids from a reservoir, etc. 

BACKGROUND OF THE INVENTION 

Patent FR-2,780,798 filed by the applicant describes a method for gradually 
deforming a stochastic model (of Gaussian type or similar) of a heterogeneous medium 
such as an underground zone, constrained by a series of parameters relative to the 
structure of the medium This method comprises drawing a number p of independent 
realizations (or representations) of the model or of at least part of the selected medium 
model from all of the possible realizations, and one or more iterative stages of gradual 
deformation of the model by carrying out one or more successive linear combinations of 
p independent initial realizations, then composite realizations successively obtained 
possibly with new draws, etc., the coefficients of this combination being such that the 
sum of their squares is 1 . 
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Patent FR-2,795,841 filed by the applicant describes another method for gradually 
deforming the representations or realizations, generated by sequential simulation, of a 
non necessarily Gaussian stochastic model of a physical quantity z in a gridded 
heterogeneous medium in order to adjust them to a series of data relative to the structure 

5 or to the state of the medium, which are collected by prior measurements and 
observations. The method essentially comprises applying an algorithm allowing gradual 
deformation of a stochastic model to a Gaussian vector with N mutually independent 
variables, which is connected to a uniform vector with N mutually independent uniform 
variables by the Gaussian distribution function so as to define realizations of the 

10 uniform vector, and using these realizations to generate representations of this physical 
quantity z, which are adjusted in relation to the data. 

The above methods are applicable to gridded models (pixel type models) suited for 
modelling continuous quantity fields and they are therefore ill-suited for modelling of 
zones crossed by fracture networks or channel systems for example. 

15 Models based on objects are spatial arrangements of a population of geometrically 

defined objects. Basically, an object type model is a Boolean model that can be defined 
as a combination of objects identical by nature with a random spatial distribution. 
Boolean models (of object type) are of great interest for the geometric description of 
heterogeneous media such as meandrous deposit systems, fracture networks, porous 

20 media on the grain size scale, vesicle media, etc. Geologic objects are defined by then- 
shape and size. Their location in the field is defined by taking account of their 
interactions : attraction-repulsion, clustering tendency, etc. 
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Unlike pixel type models, the models based on objects can provide for example 
realistic geologic representations of an underground reservoir at an early stage where 
the data obtained by in-situ measurement are still rare. 

The prior art in the sphere of object type models is notably described in the 
5 following publications : 

- Matheron, G., 1967, « Element pour une theorie des milieux poreux», Masson, 
Paris ; 

- Matheron, G., 1 975, « Random sets and integral geometry », Wiley, New York ; 

- Serra, J., 1982, « Image analysis and mathematical geology », Vol.1, Academic 
10 Press, London ; 

- Stoyan, D.S. et al., 1995, « Stochastic geometry and its applications », 2 nd Edition, 
Wiley, Chichester ; 

- Lantuejoul, C, 1997, Iterative algorithms for conditional simulations, in Baafi and 
others, eds. ; Geostatistics Wollongong 96, Vol.I, Kluwer Acad. Publ., Dordrecht, The 

15 Netherlands, p.27-40. 

The position of the objects in an object type model is distributed according to a 
Poisson point process well-known to specialists. The shape and the size of the objects 
are independent of their positions. This model can be generalized by a combination of 
objects of different nature and/or using a non-stationary density Poisson point process. 

20 Although Boolean models have been widely studied in the literature, there is no 

coherent and efficient method for constraining these models to the physical data, 
notably hydrodynamic data, which is however a major challenge for their application to 



reservoir engineering. The methods allowing gradual deformation of pixel type 
stochastic reservoir models such as those described for example in the two 
aforementioned patents cannot be directly used for Boolean models. Coristraining 
Boolean models to hydrodynamic data for example requires development of coherent 
algorithms for deformation and displacement of the objects. 

SUMMARY OF THE INVENTION 

The method according to the invention allows to generalize the gradual deformation 
technique described in the two aforementioned patents to stationary or non-stationary 
Boolean models, with or without geometric constraint as regards the well. The method 
appears to be particularly useful notably for reservoir engineers anxious to coherently 
and efficiently adjust object type reservoir models. 

The method according to the invention allows to gradually deform an initial 
realization defining the distribution of a set of objects in a zone of a heterogeneous 
medium such as a geologic structure, formed by simulation of an object type stochastic 
model, the objects being distributed in the zone according to a Poisson point process in 
form of figurative points with a point density X(x) that varies according to their position 
(x) in the zone. It essentially comprises the following two stages : 

- generating a realization of a uniform random vector according to which the position 
of each object is defined while respecting density A,(x), and 

- gradually modifying the uniform random vector according to a gradual deformation 
procedure, so as to obtain gradual migration of each object and consequently gradual 
change in the distribution of the objects in the zone, until a final realization best 



adjusted to parameters relative to the structure of the medium, such as hydrodynamic 
parameters, is obtained. 

It is possible to limit the migration of the figurative points in a subdomain of the 
zone (a well crossing the zone for example) by imposing a zero point density in the 
complementary part of the subdomain. 

According to an implementation mode, one gradually changes from a realization 
containing a first set of Nj points to a realization containing a second set of N 2 points by 
constructing a chain N(t) of Poisson numbers between the two numbers N, and N 2 , 
using the gradual deformation procedure. 

It is also possible to gradually deform the size, the shape and the orientation of an 
object during its migration by using the gradual deformation procedure. 

In cases where it contains errors, it is also possible to gradually adjust the point 
density X(x) by using the gradual deformation procedure. 

The final realization obtained by means of the method is a realistic representation of 
the configuration of the objects in the heterogeneous medium. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other features and advantages of the method according to the invention will be 
clear from reading the detailed description hereafter, with reference to the 
accompanying drawings wherein : 

- Figures 1A, IB respectively show a density function of a non-stationary Poisson 
point process and a realization of the non-stationary Poisson point process generated by 
the sequential method ; 



- Figures 2A to 2D show various examples of gradual migration trajectories between 
two points of a non-stationary Poisson point process ; 

- Figures 3A to 3H show various successive stages of a realization of the non- 
stationary Poisson point process during the gradual migration of the points ; 

- Figure 4 shows the possible domains of migration of a disk in the case of three 
conditioning points ; 

- Figure 5 shows, by way of example, a complete chain of successive realizations of a 
Boolean simulation of elliptical objects where deformation parameter t ranges between 
-7i and Ti by intervals At=0.l7i ; and 

- Figure 6 shows, by way of example, an incomplete chain of successive realizations 
of a Boolean simulation of elliptical objects where deformation parameter t ranges 
between 0 and 0.27t by intervals At=0.0l7t. 

DETAILED DESCRIPTION 

General points 

The geologic objects to which the method applies are for example more or less large 
fractures inside a reservoir zone, or channels. The method can also apply to granular or 
vesicle type structures of much smaller size. All these objects are difficult to model by 
means of pixel type models. 

The progressive migration operations described hereafter apply to an initial model 
where the positions of the objects are represented by point configurations (referred to as 
point procedures). The distribution of these points varies in density according to then- 
positions in the modelled zone. This distribution is based on various data known from 



measurements or observations : geomechanical measurements obtained in wells for 
example, seismic data obtained from prior seismic operations. 

Starting from this initial model, the method allows to gradually deform the initial 
distribution according to a series of rules so that, in the final distribution, the model is 
optimized in order to better correspond to constraints imposed on one or more physical 
parameters such as, for example, a distribution of permeability values. The migration 
rules are such that a great number of different points of the model can be globally 
displaced from a limited number of control parameters. 

An objective function measuring the difference between the physical parameters 
from the real heterogeneous medium and those simulated on a realization of the 
stochastic model is generally defined. The value of the objective function therefore 
depends on control parameters of the stochastic model. The values of these control 
parameters are obtained by minimizing the objective function. 

Summary of the Poisson point process 

The Poisson point process is a countable random set of points distributed in any 
space 9? n . This set of points has the following characteristics : 

Let D be a domain of 9l n . If the volume of D, denoted by I D I, is finite, then the 
numbers of points lying in D, denoted by N(D), follows a Poisson' s law of parameter 
X ID I. Let : 

? [ #(Z )) = ^lM Vn >o (i) 



where X is called the point process density ; it measures the average number of points 
lying in a unitary volume domain of 9T. 

Let D b D 2 ,..., D k be domains of 9t a unconnected two by two, then the numbers of 
points lying in these domains NODj), N(D 2 ),..., N(D k ) are mutually independent random 
variables. 

Conditionally on N(D)=rip, these rip points are independent and evenly distributed in 

D. 

The description hereafter presents the Poisson point process in the finite domain D. 
Migration of a stationary Poisson point process 

Consider the problem of the migration of a realization of a stationary Poisson point 
process in D rectangular. In order to simplify the presentation, D is assumed to be a 
unitary hypercube [0,1 f with n dimensions. Let X! and x 2 be two independent points 
uniformly drawn from [0,1 ] n . We define a trajectory between x x and x 2 by : 

x(t) = g[g _1 (Xj)cos? + G' 1 (x 2 )sin t] (2) 

where G is the centered and reduced Gaussian distribution function. According to the 
gradual deformation algorithm described in the aforementioned patent FR-2,780,798, 
for any t, x(t) is a uniform point in [0,1]°. When the two points are fixed, the trajectory 
of the gradual migration between them is entirely determined. The place change of one 
of the two points will change the migration trajectory. We can show that the trajectory 
defined by Equation (2) is symmetrical in relation to the centre of the domain [0,1 ] n , 
whatever the number of dimensions n. This suggests that, even if the two points are 



isolated in a corner/side of the domain, the trajectory of the gradual migration between 
them can always reach the opposite part of the domain. 

Migration of a non-stationary Poisson point process 

Consider now the Poisson point process in a domain D of density A,(x) of general 
shape. The number of points in D is a random Poisson variable of average 
2(D) = j b JL(x)dx . These points are independently distributed in D according to the 
probability density function : 
fix) = Z(x)/A(D) xeD (3) 

Simulation of a Poisson point process of density X(x) in D can be carried out in two 
stages : 

- generating a number n according to Poisson' s law of average equal to X(D), then 

- generating n points in D irrespective of one another according to the same 
probability density f(x). 

The example of Fig.lA shows a density function, and Fig.lB shows a realization of 
a Poisson point process of density A,(x). 

If law f(x) is simulated by inversion of its distribution function, then point x 
corresponds to a uniform vector u. The gradual migration algorithm can thus be applied 
to the Poisson point process of density A.(x). By construction, this method preserves the 
density and the number of points of the initial process. Figs. 2A to 2D show four 
migration trajectory examples. 
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Take the example of the non-stationary Poisson point process in a two-dimensional 
domain D = [0,1] 2 whose density increases linearly in axis X and is constant in axis Y. 
Let (x,y) be the vector of the coordinates of a point in D of this process. Then, (x,y) 
admits the bivariable probability density : 

f(x,y) = 2x, (x,y)e[0,if (4) 

5 

Simulation of a point according to the above law (4) is simple : 

- generating abscissa x according to the linear density f(x) = 2x, then 

- generating ordinate y uniformly between 0 and 1 . 

Let : 

x = 4u 

JP = v (5) 
10 where u and v are two independent and uniform numbers between 0 and 1 . 

We can thus apply gradual migration algorithm (2) to the uniform vector (u,v) in 
order to establish a migration trajectory for point (x,y) in D. 

In practice, density function X(x) often comes in form of a grid. By way of example, 
consider the case of a two-dimensional grid of M x N nodes. Let Xj and y- s be the 
1 5 coordinates of node (Lj). The marginal law for x ; is then : 

/(*,)= ^ 

and the conditional law for y j5 knowing x ; , is : 

( 7 ) 
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The points of the Poisson process of density l(x) can thus be readily generated on a 
grid from these mono variable probabilty laws. 

By way of illustration, we have constructed a two-dimensional density function 
from a centred and reduced Gaussian simulation. The variogram is anisotropic and of 
5 Gaussian type. The principal anisotropy directions are diagonal in relation to the 
coordinates system. The scale factors in the principal anisotropy directions are 0.3 and 
0.1 respectively. The size of the field is lxl and it is discretized in 1000x1000 pixels. 
The Gaussian numbers are changed into positive numbers according to the following 
expression : 

l(x) = 400e^ x) (8) 

10 where Y(x) is the Gaussian simulation. Fig.lA shows the density function thus 
constructed. Fig. IB shows a Poisson point process realization admitting the density 
function of Fig.lA. 2000 points are generated according to the sequential method. 

Figs. 3 A to 3H illustrate the evolution of a Poisson process realization during 
gradual migration of the points. It can be seen that the point density is respected during 
15 the migration. 

Migration of conditional Boolean simulation ob jects 

The above method can be immediately applied to the migration of a domain S of 
any form in D. In fact, the migration of a point in S can be carried out in D using the 
truncated probability density function : 

f s (x) = A,(x)l xeS /A,(S% xeD, SaD (9) 
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Using the Markovian iteration algorithm described in the reference to Lantuejoul, 
1997, it is possible to simulate a Boolean model in a domain D knowing that two 
subsets Ci and C 0 of D respectively belong to the union of the objects and to its 
complement. Then, the migration algorithms in any domain can be used for gradual 
5 deformation of the Boolean simulations conditioned by geometric well data. In fact, 
from a conditional Boolean simulation and without compromising conditioning by Q 
and C 0 , the objects must move only in their respective domains defined according to the 
shape of the objects and the configuration of d and C 0 . 

Consider an object A of a conditional Boolean realization, which includes a subset 
10 Cia of d and excludes C 0 . If, during migration of A, it still has to include C 1A but 
excludes C 0 , the allowed migration domain for object A is : 

D A = {x:C 0 nA x = 0;C U a A x ] (10) 

Fig.4 shows a case with three conditioning points and the eight possible domains of 
migration of a disk. If, for example, the disk is supposed to always cover points (a) and 
(b) but avoid point (c), then its centre can move only in domain 3. 

15 Migration with appearance and disappearance of points 

The number of points in D of a point process follows a Poisson's law of parameter 
1(D). It is therefore necessary to vary the number of points in D during their migration. 
In this section, we first present two methods of constructing Poisson number chains, 
then we shall see how to migrate between two sets of points whose cardinals are not 

20 identical. 
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We try to construct a chain of Poisson numbers between two Poisson numbers Ni 
and N 2 , independently generated by inversion of the distribution function. Let Ui and U 2 
be two independent uniform numbers (between 0 and 1) from which numbers N, and N 2 
are obtained. 

N 2 = F-\U 2 ) (11) 

where F 1 represents the inverse distribution function of Poisson's law. According to the 
gradual deformation algorithm, we can construct a chain of uniform numbers between 
Ui and U 2 by : 

UiO^GlG-'OJOcost + G-'iUJsint] (12) 
By inversion of the distribution function of Poisson's law, we obtain a chain of 
Poisson numbers. 

N(t)= F~ l [U(t)] (13) 
The inverse distribution function can be calculated by dichotomy. However, if the 
parameter of Poisson's law is too great, this method remains quite costly. 

In order to save calculating the inverse distribution function of Poisson's law, 
another way of generating Poisson's numbers can be envisaged. It is well-known that 
the number of jumps of a Poisson process of parameter 1, in a segment of length X, 
precisely follows a Poisson's law of parameter X. Simulation of a Poisson process is 
easy to do. In fact, the intervals between two consecutive jumps of the process are 
independent and follow the exponential law of parameter 1. By gradually deforming the 
exponential segments, we obtain a chain of Poisson numbers. It is easy to construct a 
chain of exponential numbers because calculation of the inverse distribution function of 
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the exponential law is easy. It can be reminded here that the distribution function of the 
exponential law of parameter 1 is written as follows : 

F(s) = l-e~\ s>0 (14) 
and its inverse function is : 

F-\u) = -\n(l-u), 0<u<l (15) 

The fact remains that, since the number of jumps in the segment of length A, is X on 
average, the number of chains of exponential numbers is close to X on average. 

Consider now the migration of a set of N x points (set 1) to a set of N 2 points (set 2). 
Since N! * N 2 , the migration from set 1 to set 2 necessarily implies the appearance or 
the disappearance of certain points. The number of points that will appear or disappear 
is determined by the chain of Poisson numbers between N x and N 2 . The migration 
algorithm is as follows : 

a) Calculate the maximum number of chain N(t). 

b) Complete set 1 with - N t points and set 2 with - N 2 points. 

c) For each set, order the points from 1 to N^. All the initial points are ordered 
from the beginning. 

d) Calculate the migration trajectory of point n of set 1 to point n of set 2 

e) For each state t of the set, remove the - N(t) last points. 



